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SUMMARY 


This report presents the research performed at the Catholic 
University of America on the research grant entitled: “Robust 
Design of Distributed Controllers for Large Flexible Space 
Structures, " under the Grant No NAG 5-949, between June 30th, 1987 
and June 30th, 1988. 

Independent Modal Space Control (IMSC) method avoids control 
spillover generated by conventional control schemes such as Coupled 
Modal Control by decoupling the large flexible space structure into 
independent subsystems of second order and controlling each mode 
independently. The IMSC implementation requires that the number of 
actuators be equal to that of modeled modes, which is in general 
very huge. Consequently the number of required actuators is 
unreal izable. 

In this report two methods are proposed for the implementation 
of IMSC with reduced number of actuators. In the first method, the 
first m modes are optimized, leaving the last (n-m) modes 
unchanged. In the second method, generalized inverse matrices are 
employed to design the feedback controller so that the control 
scheme is subopt tmal with respect to IMSC. The performance of the 
proposed methods is tested by performing computer simulation on a 
simply supported beam. Simulation results will be presented and 
discussed. 


V 1 ' ~ A 

Ft* > • 

\<Yv- *■ ■■■/• 

‘ r; ro i\uc \ v* *> a V '■ - 


** c-t vj ^ n c vt - 

C' \ rv- '\ ^ ' 

f ^ ~ f * r. ► "Av ' 

r & k\r *0 > '■ ■ « 


t INTRODUCTION 

The development of the space shuttle has opened the 
possibility of constructing very large space structures (LSS) In 
space for space explorations. Two control problems for LSS are 
attitude control and shape control. Complex missions Impose many 
stringent requirements on shape and attitude of the LSS, which lead 
the control researchers to the concept of distributed active 
control that places on the structure a number of sensor /actuator 
pairs in oder to optimize the LSS performance and behavior. Active 
control of LSS has been an active research area in the last several 
years 11 3— [73. A large number of control schemes has been 
developed for LSS, but they represent one form or another of modal 
control [1]. Two main modal control schemes are the Coupled Modal 
Control (CMC) and the Independent Modal Space Control (IMSC). The 
CMC employs an active controller that consists of a state estimator 
and a state feedback while the IMSC decouples the LSS into n 
Independent subsystems according to n controlled modes and controls 
each mode Independently by means of a modal filter [5] and an 
optimal controller. It is well-known that CMC causes control and 
observation spillover, which together can destabilize the LSS [1]. 
IMSC eliminates control and observation spillover because each mode 
is controlled independently. However the implementation of IMSC 
requires that the actuator number be equal to the number of modeled 
modes, which is usually very huge for the modeling of the LSS to be 
faithful. This fact presents a fundamental limitation of IMSC for 
the required number of actuators is impractical. The main 
objective of this report is to propose methods implementing the 
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IMSC with a milder requirement of the actuator number. In 
particular, we wl 1 1 develop two control schemes that uses a reduced 
number of actuators to control all modeled modes In such a way that 
the closed- loop system modes are as Identical as possible to the 
optimal modes specified by the IMSC scheme. In the first control 
scheme, the first m modes are optimized, leaving the last (n-m) 
modes unchanged. In the second scheme, generalized Inverse 

matrices are employed to design the feedback controller so that the 
control scheme is subopt Imal with respect to IMSC. 

Matrix notations used In this report are given below: 

i 0 0 

M . . . 0 
2 

0 M 

0 = mxn null matrix 
mxn 

1 = nxn identity matrix 
n 

Z SUMMARY OF INDEPENDENT MODAL SPACE CONTROL 

A large flexible space structure can be described by the 
following partial differential equations [4]: 

M(P) 3 2 u(P,t)/9t 2 + Lu(P, t) = f(P,t) (1) 

that must be satisfied at every point P of the domain D, where 
u(P,t) Is the displacement of point P, L a linear differential 
self-adjoint operator of order 2p, expressing the system stiffness, 
M(P) the distributed mass, and F(P,t) the distributed control 
force. The displacement u(P,t) Is subject to the boundary 
conditions: 

T.u(P,t) =0; 1 = 1,2 p (2) 
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where T , 1 = 1,2 p are linear differential operators of order 

ranging from 0 to (2p-1). 

The associated eigenvalue problem Is formulated by: 

L® (P) = A M(P) * (P); r = 1,2,... (3) 

r r r 

with the boundary conditions: 

T i * r (P) =0; 1 = 1,2 p; r = 1,2,... (4) 

where A^ Is the rth eigenvalue and 4^(P) Is the eigenfunction 
(sometimes also known as Mode Shape) associated with A^. Suppose 
the operator L is self-adjoint and positive definite, and all 
eigenvalues are positive and are ordered so that A^<A <. . . . Since 
L Is self-adjoint, the eigenfunctions are orthogonal and therefore 
can be normalized such that: 


fM$ t dD = 6 
D J r 8 

and f $ L* dD = A 8 ; r,s = 1,2, . . 

rJ * r r rs 


(5) 

( 6 ) 


where 6 Is the Kronecker Delta, 
rs 

Using the expression theorem [3], the solution of u(P,t) can 

be obtained as: 

00 

u(P,t) = Z * (P)u (t) (7) 

. r r 

r»1 

where u (t) Is the modal coordinate. Substituting (7) Into (1), 
r 

multiplying both sides of the resulting expression by 9 , 

s 

integrating over D and employing (5) and (6), we obtain 

u (t) + « 2 u (t)=f(t); r = 1,2, ... (8) 

r r r r 

In (8), the mode (or natural frequency) is defined as 

u = (A ) 1/2 ; r = 1.2 (9) 

r r 

and the modal control force f^(t) is computed by: 
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( 10 ) 


f (t) = f * (P)f(P,t) dD 

r h r 

In practice, the Infinite series In (7) Is truncated as 
n 

u(P,t) = I * (P) u (t) (11) 

r«1 r r 

where n Is chosen to be sufficiently large so that u(P,t) can be 
represented with good fidelity. In this case we are dealing only 
with the first n modes. 


Eq. (8) 

can be 

transf ormed 1 nto 

state equation form 

fol lows: 

x(t) • 

= A x(t) + 

W(t) 

(12) 

where x(t) = 

(xj(t) x 

l( t) ....x T (t)] T 
2 n 

(13) 

W(t) • 

= (wj(t) 

w*(t) w T (t)] T 

2 n 

(14) 

A = Block d1ag(A ,A , A^) 

(15) 

x (t) 

= [u (t) 

u (t)/o ]* T 

(16) 

r 

r 

r r 


W (t) 

= [0 f (t)/o ] T 

(17) 

r 

r 

r 


and A = 

r 

0 u 

r 

-o 0 

l r J 


(18) 


for 1=1,2 n. 


REMARK 1: 

At this point we should distinguish between the eigenvalues of 

A and the eigenvalues of the flexible space structure. From (15) 

and (18) the 2n eigenvalues of the matrix A that represent also the 

open-loop poles of the system consist of n pairs of imaginary 

numbers: +oi, +uj, ..,+«•> j. However according to (9), the 

I 2 n 

2 2 

eigenvalues of the flexible space structure are given by iy ty 

o 2 . We also note that the rth mode u is the magnitude of the 
n r 
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eigenvalues of A^. In other words, a mode corresponds to 2 
eigenvalues of one subsystem of the complete flexible system. 

The system described by (12) consists of n subsystems given 
by 


x (t) = A x (t) + W (t); r = 1,2 n (19) 

r r r r 

The essence of IMSC is to choose VMt) such that It depends 
on x (t) alone. Thus 


W r (t) = G x ; r = 1,2, 
where C = K" 

r L9 rJ , 9,J 

are (2x2) gain matrices. 


. ,n 

r = 1,2 n 


(20) 

( 21 ) 


Substituting (16) and (17) Into (20), we find that G must 

r 


assume the following form: 


r o o' 
r Lg r21 g r 22 J 

For optimal control, g .. and g „ 

r 2 1 r 22 

the following quadratic cost function 
problem): 


1,2 n (22) 

should be determined such that 
is minimized (linear regulator 


J 


n 


= Z J 

r=1 


r 


(23) 


where J = f (x T Q x + W T RW )dt 

rj rrr rrr 

0 

and are positive semi definite and positive 
matrix, respectively, associated with the rth mode. 

The form of G^ given by (22) requires that 
given below [7]: 


(24) 

definite weighting 

R assume the form 
r 




r = 1,2 n. 


(25) 


Since depends on x^ alone as seen in (20), J can be minimized by 


S 


minimizing each J , Independently. From optimal control theory 


[15], the optimal solution for is given by 
G (t) = R 1 K (t); r = 1,2 n 


(26) 


where K (t) Is the solution of Rlccatl equation: 

K (t) = -K A - A T K + KR* 1 -Q; r = 1,2 n (27) 

p rr rr rr r 

with boundary condition K^(T) = 0. 

Selecting Q = o 2 I and substituting (25) and (18) into (27) 

r r 2 

yield 


K, = 2o K + r' 1 K 2 - o 2 

11 r 12 r 12 r 

+ r_1 < 

12 r 22 11 r 12 22 

* -12 2 
K = -2o K + r K ~ u 
22 r 12 p 22 r 


(28) 

(29) 

(30) 


where K (1, j = 1,2) are the elements of K^, and = K^.The 


steady state solutions of (28) -(30) can be found by letting K 


11 


K 12 = * 22 = °- 


* 2 

Introducing R = (l/o )R where 
r r r 


* 0 ] 

'r - 0 r* j 

L r J 


r = 1,2 n 


the steady-state solutions of (28) -(30) are obtained as 

„ t, 2 , * 1/2. 

K = K = o ( -o r + a ) 

12 21 p p p 


22 


/[r - 2u 2 r* 2 + 2o r* a 1/2 J ,/2 

r r r r r r 


K ■ o 2 [1/w 2 + (2/or*)a 3/Z - 2r*r 2 -r* ] 1/2 
11 r r r r r r 

2 *2 * 

where a = u r + r . 

r r r 

Substituting ( 32 ) - ( 34 ) Into (26) yields 


(31) 

ft 

(32) 

(33) 

(34) 

(35) 


G = 

p 


o - b -[2o (-o + b) + r 

r p r r 


(36) 


where b * (o 2 + r** 1 )^ 2 
r p 

From (17), (20), (22) and (35), we obtain 


(37) 


f (t) ■ (o - b)u u (t) - [2 o (-u + b) + r*' 1 ] 17 ^ (t) (38) 

p ppp rp pr 
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for r * 1,2 n. 

Substituting (38) into (8) yields 

u (t) + [2u (b - o ) + p*" 1 ] 1/2 u (t) + o bu (t) = 0 (39) 

r r r r r r r 

From (39) the closed- loop poles (eigenvalues) of the system are 
S = =-1/2 [2u (b-o ) + r** 1 ] 1/2 + 1/2 [-2o (b+w ) + r*'V /2 (40) 

r1,r2 r p r p r p 

for r = 1,2 n. 

The closed- loop system is stabilized since the poles given in 
(40) are either complex conjugates with negative real parts, or 
both real and negative. Figure 1 illustrates the concept of 
Independent Modal Space Control in the state equation form. As we 
observe, each mode is controlled Independently by an actuator. 
Therefore the number of required actuators must be equal to the 
number of modes that is in general very high. This facts 
represents a principal disadvantage of IMSC. On the other hand, 
one main advantage of IMSC is the elimination of control and 
observat Ion spill over . 

3. IMPLEMENTATION OF INDEPENDENT MODAL SPACE CONTROL 

The IMSC scheme is represented by Equations (19) and (20) that 
however are In the modal form. In order to Implement the IMSC, a 
finite number of discrete sensors and discrete actuators must be 
employed. 

DISCRETE SENSORS AND NODAL FILTER 

Suppose there are 2p discrete sensors consisting of p 
displacement sensors and p velocity sensors located at p locations 
Using (11), the output of the displacement sensor and 

12 p 
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velocity sensor can be expressed as 


y. = u(P.,t) = I * (P.) u (t), 1 = 1,2 p 

i i . r i r 

r*1 
n 

and y, = u(P.t) = E # (P.) u (t), 1 = 1,2 p 

f i . r i r 

r= 1 


( 41 ) 


y(t) = [y, y, y, y, •••y„y„]'. 


If we define an output vector y(t) containing both 
displacement and velocity sensor output as 

(42) 

n 

then using (41), we can write 

y(t ) = C x ( t ) (43) 

where C = [c ], 1 , J = 1,2 p (44) 


and 


ii 


*.(P.) 

j » 


o.*.(P.) 

J J » J 


(45) 


In order to Implement (20) modal state vector x^(t) must be 
extracted from the sensor outputs. A device called Modal Filter 
[5] Interpolates the discrete sensor measurements to obtain 

A X 

continuous displacement and velocity profiles u(P,t) and u(P,t) and 

then computes the estimated modal displacement and velocity using 

the expansion theorem stated by (7): 

u (t) = l*M(P) * (P) u(P,t)dD (46) 

r r 

u (t) = f*M(P) * (P) u(P,t)dD (47) 

r J r 


for r = 1,2, . . ,n 


DISCRETE ACTUATORS 

Since it Is Impossible to control force at every point In the 
domain D, the distributed control force is realized by m (m<n) 
discrete point force actuators applied at m points P , P ? P^ 
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( 48 ) 


In the domain D as given below 

■ 

f(P,t) = z SCP-P^ F.(t) 

i * 1 

where S(P-P.) Is a spatial Dirac Delta function and F.(t) Is the 

force applied by the 1th actuator on the point P . 

Now substituting (48) Into (10) yields 

f (t) = I* (P) Z «(P-P.) F.(t) dD (49) 

r Jr ii 

D 1*1 

From the property of the Dirac Delta function, (49) can be reduced 
to 

* 

f (t) * z * (P.) F (t) (50) 

p . „ r 1 i 

i * 1 

If we define a force vector F(t) such that 

F(t) = [ F^t) F 2 (t)....F B (t)] T , (51) 

then using (50) the relation between F(t) and W(t) can be expressed 
by 


W(t) = B F(t) 


(52) 


where 


B = [B T B T ....B T ] T 
12 2n 

W°W 1 " 1 ’ 2 " 

*, (P 2 )/U 1- 

for 1 = 1,2 n. 


*,(P )/o 1 

i it 1 


(53) 

(54) 

(55) 


4. PROBLEM STATEMENT 

Since In this paper, the Implementation problem Is focused on 
actuators, we suppose that there exists a modal filter which takes 
2p measurements of p displacement sensors and p velocity sensors 
and produces a "perfect" estimated state x(t) of x(t). It Is 
obvious that If the modal displacement u^(t) and modal velocity 
iMt) can be perfectly estimated as shown by (46) and (47), then In 
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view of (16) and (13), the state x(t) can also be perfectly 
estimated. Here the state x(t) Is said to be perfectly estimated 
If x(t) approaches x(t) as t — > «. In control theory, the 
corresponding state estimator Is said to be asymptotically stable 
[ 10 ]. 

Figure 2 Illustrates the Imlementatlon of IMSC in which the 
optimal state feedback law is denoted by 



W(t) = G x ( t ) 



(56) 

where 

W(t) = [W^t) 

w 2 T (t).. 

..W T (t)] T 
n 

(57) 

and 

G = Block diag 

(( V V 

...G n ) 

(58) 

In 

section 3, the 

opt 1 ma 1 

solution f or 

IMSC Is obtained 


modal space, namely from (20) 

W(t) = G x ( t ) (59) 

Under the assumption of perfect state estimation, we note that 

A 

W — > Gx(t) as t-» oo. Therefore the optimal solution can be 
achieved If (59) Is satisfied. To make W(t) equal to Gx(t), the 
matrix D In Fig. 2 Is chosen such that W(t) = W(t). 

From Fig. 2 we note that 

F(t) * D W(t) (60) 

Substituting (60) Into (52) yields 

W(t) = B D W(t) (61) 

Obviously to make W(t) = W(t), D Is chosen such that 

00*1, (62) 

2n 

From the structure of B as given by ( S3 ) — ( 55 ) , each (21—1 ) th 
row (odd row) of BO, for 1=1, 2,..., n Is a row of zeros. We realize 
that (62) can never be satisfied. However, noting that each odd row 
of W(t) Is also a row of zeros, If we define 
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( 63 ) 


B(t) = [B* B* ..B* ] T 
2 4 2n 

and D(t) = [ d\ D ] D* ] T (64) 

2 4 2n 

where B. Is the 1th row of B and the 1th column of D, then 

choosing a matrix D such that 

B D = I (65) 

n 

will ensure that W(t) = W(t). It Is noted that If (65) holds, then 
BO Is a modified Identity matrix of order (2nx2n) whose main 
diagonal elements are 0 at the (21-1,21-1) position and 1 at the 

(21, 2i) position for 1 * 1,2 n. 

The solution of (65) for D needs the Inverse of B. However, 
the existence of the Inverse of B require that the number of 
actuators be equal to the number of modeled modes (m=n). This fact 
presents the principal limitation of IMSC because the number of 
modeled modes Is usually very high, resulting In an unrealizable 
number of actuators. 

The problem considered in this report is formulated as to 
develop a control scheme that allows one to use a reduced 
number of actuators to control all modeled modes such that the 
closed- loop system is as optimal as possible compared to IMSC. 

In the following we will propose two control schemes that will 
accomplish the above objective. 

5. THE FIRST CONTROL SCHEME 

The development of the first control scheme Is represented In 
the following theorem: 
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THEOREM 1: 


Consider a large flexible space structure whose description 
and solution are given by (1) and (11), respectively. If the 
operator L Is self-adjoint and the state x(t) Is perfectly 
estimated, then there exists a control scheme with m actuators, 
which optimizes the first m modes with respect to the cost function 
(24) and leaves the last (m-n) natural modes unchanged, thus 
ensuring the system stability. 


Proof: 

If the hypothesis of Theorem 1 Is satisfied, then employing m 

actuators placed at m distinct points P , P P in the domain 

12 m 

D results In a matrix B that can be partitioned as 
B. 


B ■ 


1 

§ 


(66) 


where B^ls an (mxm) matrix and B 2 an (n-m)xm matrix. Based on the 
property of the rows of B specified by (55), without loss of 
generality we can assume that the first m rows of B are linearly 
Independent for the locations of m actuators are distinct. Thus § 
Is a nonsingular matrix. Now If a matrix D Is chosen such that 

] (67) 


6 = [B' 1 0 

1 mx(n-m) 


then It Is clear that 


BO 


I 


s 2 § ;’ 


pnx(rvm) 


(68) 


(n-m)x(n-m) J 

We proceed by considering the closed- loop feedback system given In 
Fig. 2 that can be described by 

x(t) - A x ( t ) + B D G x ( t ) (69) 

under the assumption of perfect state estimation. In the steady 
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state (t-» oo), (69) can be rewritten as: 

x(t) = (A + BOG) x ( t ) (70) 

Now If 5 Is chosen as in (67), then 

A+BDG = Block d1ag(A ,A A ) + 

12 n 

I o 

+ 2* 2mx2(n-*) 

X 0 

2(n-m)x2(n-*)- 

where X Is an 2(n-m)x2m matrix whose elements are obtained from 
— — * 

BB and I is a modified (2mx2m) identity matrix whose main 
2 12 * 

diagonal elements are 0 at the (21-1,21-1) position and 1 at the 

(21,21) position for 1 = 1,2 m. 

Recalling the form of G for r = 1,2 n as given In (22), 

r 

we can rewrite (71) as 


Block diag(G.,G„, G) (71) 

12 n 


A + BDG = 


Block dlag (A + G , A + G , ...,A+ G ) 

i H 2 2 fl® w 


2m x2 ( n - m ) 

Block diag(A , A , ...» A ) 
12 n 

(72) 


From (72) the closed- loop eigenvalues can be computed as 
m ■ 

<r(A + BDG) = U <r(A. + G) + U <r( A. ) (73) 

i.i ■ 1 ..i 1 

* 

where U S. denotes the union operation on the sets S., 1=1, 2,..., m. 
i «1 

and <r(A) denotes the set of eigenvalues of the matrix A. 

According to Remark 1. It is obvious from (73) that the first 
m modes are optimized and the last (n-m) modes are unchanged. 
Since optimal modes correspond to stable closed- loop poles as 
stated by (40) and natural modes correspond to pairs of imaginary 
complex conjugate eigenvalues, as pointed out In Remark 1, we 
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conclude that the above control scheme also ensures the system 
stability. We therefore just complete the proof of Theorem 1. 

REMARK 2: 

If the control scheme proposed In Theorem 1 Is Implemented, 
then there are some unwanted excitations of the (n-m) modes, which 
Is well-known as control spillover [1]. However the oscillations 
caused by exciting the (n-m) modes are Insignificant because the 
mode amplitudes tend to decrease as the number of modes is 
Increased because higher modes require more energy to excite. We 
also note that in (1), we assume the worst case where no natural 
damping exists. In practice, existing natural damping In the 
structure can suppress the unwanted oscillations of the higher 
modes. 

6. THE SECOND CONTROL SCHEME 

In this section, another control scheme Is developed using the 
concept of generalized Inverse matrices. We first present the 
following lemma: 

LEMMA 1: 

Consider the following equation: 

W(t) = B F(t) (74) 

where W(t) and F(t) are matrices consisting of even rows of W(t) 
and F(t), respectively. If the (nxm) matrix has rank m, then the 
solution for (74), which minimizes the weighted norm of error 
| |e(t) | || - | |W - BF| || * (W - BF) T S(W - BF) (75) 
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Is given by 

F(t) = B* W(t) (76) 

where the matrix 

B* = (B T SB)‘ 1 B T S (77) 

is the generalized Inverse of B . A proof of Lemma 1 can be found 
In 1141. 

The following theorem will present the development of the second 
control scheme: 

THEOREM 2: 

Consider a large flexible space structure whose description 
and solution are given by (1) and (11), respectively. If the 
operator L Is self-adjoint, then there exists a control scheme with 
m actuators (m < n) that Is suboptlmal with respect to (24) In the 
sense that the closed- loop eigenvalues are assigned as closed as 
possible to those optimal eigenvalues specified by IMSC. 

Proofs 

A control scheme with a reduced number of actuators would be 

optimal If 6 could be selected to be a right Inverse of B In (65). 

Here B Is assumed to have rank m since the discrete actuators apply 

point forces at m distinct points P , P P . From [14], It Is 

well-known that an (nxm) matrix (m < n) having rank m does not have 

any right Inverse. According to Lemma 1, because F(t) as given In 

_ * 

(76) minimizes (75), selecting a matrix D = B will minimize the 

difference §6 - I , making (65) satisfied as well as possible. 

n 

_ * 

Therefore selecting D ■ B also makes the closed- loop eigenvalues 
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as Identical as possible to those specified by IMSC. Thus there 
exists a control scheme with a reduced number of actuators, which 
Is suboptlmal with respect to (24). We just complete the proof of 
Theorem 2. 

7. COMPUTER SIMULATION STUDY 

In order to evaluate the performance of the proposed control 
schemes, we consider the control of a simply-supported beam whose 
dynamics Is given by the Euler-Bernoul 1 1 partial differential 

equation: 

g2 

El ---- u(x,t) + m ---- u(x, t) = f (x,t) (78) 

3x* 3t Z 

where for simplicity we set the mass m, the moment of Inertia I, 
the modulus of elasticity E and the length of the beam to unity. 


The boundary of the simply-supported beam are 

u(0, t) = u(1,t) = 0 (79) 

-2 -2 

---- u(0, t) = ---- u(1,t) = 0. (80) 

3x 3x Z 

The solutions for the eigenvalue problem are given by: 

X k = (kit) 2 (81) 

= (2) 1/Z s1n (kirx) (82) 


Computer simulation was performed to evaluate the second 
control scheme and to compare Its performance to that of IMSC. 
Suppose we consider 20 modeled modes and divide the whole length of 
the beam Into 20 sections, specified by x(k) = k/21 for k = 

1,2 21. Then starting with one actuator we performed computer 

simulation for different locations of the actuator. After that 
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computer simulation was performed with Increasing number of 
actuators. The simulation results are summarized and discussed 
below: 

(a) The Number of Actuators: 

Simulation results show that the dynamics of the closed- loop 
system Is affected considerably by the number of the actuators. 
As expected, the number of stable closed-loop eigenvalues Is 
proportional to that of the actuators. The maximum number of 
closed- loop eigenvalues when using 1, 2, 3, and 4 actuators is 5, 
18, 19, and 20, respectively, as showed from computer simulation. 

(b) The Location of Actuators-. 

Results show that the actuator location also affects the 
closed-loop system dynamics. Figure 2a illustrates the 
relationship between actuator location and the number of stable 
closed- loop eigenvalues for the case of 1 actuator. We note that 
the optimal actuator locations center around both ends of the beam. 
The number of the closed- loop eigenvalues decreases as the actuator 
moves toward the center of the beam. The above observations are 
consistent with the cases of 2, 3, and 4 actuators. 

(c) The Second Control Scheme versus IHSC-. 

The open- loop eigenvalues of the LSS , the closed- loop 
eigenvalues of IMSC implemented with 20 actuators and 2 actuators 
and the second control scheme Implemented with 2 actuators are 
presented In Table 1. The table shows that the second control 
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scheme Implemented with a reduced number of actuators assigns 
eigenvalues that are very close to those specified by the optimal 
IMSC. 

Figures 3 and 5 represent the vertical displacements of the 
beam center when being excited by a unit Impulse for IMSC with 20 
actuators and the second control scheme with 2 actuators, 
respectively while Figure 4 represents the case of IMSC with 2 
actuators. We observe that with a reduced number of actuators (2 
actuators), the second control scheme provides with the same 
performance as that of IMSC when implemented by 20 actuators. On 
the other hand, as we see In Figure 4, when being Implemented with 
a reduced number of actuators (2 actuators), the IMSC suffers from 
some oscillations. Beam movements for the case of IMSC with 20 
actuators and 2 actuators are Illustrated in Figures 6 and 7, 
respectively. The case of the second control scheme with 2 
actuators Is presented In Figure 8. Similar to the previous case, 
the performance of the second control scheme is better than that of 
IMSC when being Implemented with a reduced number of actuators. 

8. CONCLUSION AND FUTURE RESEARCH 

In this report, we first reviewed the theory of Independent 
Modal Space Control In the context of optimal control. The 
requirement that the number of actuators be equal to that of 
modeled modes for the Implementation of IMSC was pointed out as a 
principal limitation. Assuming that the large space structure Is 
self-adjoint, we developed two control schemes Implemented with a 
reduced number of actuators, which are suboptlmal with respect to 
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the IMSC performance. The first control scheme optimizes the first 
m modes leaving the last (n-m) modes uncontrolled. The second 
control scheme employs the method of generalized Inverse matrices 
to assign closed-loop eigenvalues, which are as close as possible 
to those specified by IMSC. Computer simulation performed on a 
simply-supported beam showed that the second proposed control 
scheme performs better than the IMSC when both are Implemented by a 
reduced number of actuators. Future research should focus the 
attention on the optimization of the number of actuators and the 
placement of actuators In order to minimize the control spillover. 
The effect of modal filter In terms of observation spillover should 
be Investigated. The method of arbitrary eigenvalue assignment [8] 
should also be considered for the design of control schemes for 
large space structures. 

The research performed under this grant were transmitted to 
one semiannual progress report [11], one conference paper [12] one 
published journal paper [9] and one technical paper to be submitted 
to a major journal [13]. 
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Figure 2a: Stable Mode Number versus Actuator Location 
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CLOSED-LOOP EIGENVALUES 


Ta bl e 1 : Open-Loop Eigenvalues and Closed-Loop Eigenvalues of 
Various Control Schemes 
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. l~ i 9 u re — 6 : Beam Response When the IMSC is implemented by 20 Actuators 
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ure 7 : Beam Response When the IMSC is implemented by 2 Actuators 







rafter 2 sec 



a tors 





